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We review advances in the modeling of protoplanetary disks. This review will focus on the 
regions of the disk beyond the dust sublimation radius, i.e. beyond 0. 1 - 1 AU, depending on the 
stellar luminosity. We will be mostly concerned with models that aim to fit spectra of the dust 
continuum or gas lines, and derive physical parameters from these fits. For optically thick disks, 
these parameters include the accretion rate through the disk onto the star, the geometry of the 
disk, the dust properties, the surface chemistry and the thermal balance of the gas. For the latter 
we are mostly concerned with the upper layers of the disk, where the gas and dust temperature 
decouple and a photoevaporative flow may originate. We also briefly discuss optically thin 
disks, focusing mainly on the gas, not the dust. The evolution of these disks is dominated by 
accretion, viscous spreading, photoevaporation, and dust settling and coagulation. The density 
and temperature structure arising from the surface layer models provide input to models of 
photoevaporation, which occurs largely in the outer disk. We discuss the consequences of 
photoevaporation on disk evolution and planet formation. 



1. INTRODUCTION 

Dusty circumstellar disks have been the focus of intense ob- 
servational interest in recent years, largely because they are 
thought to be the birthplaces of planetary systems. These 
observational efforts have yielded many new insights on 
the structure and evolution of these disks. In spite of ma- 
jor developments in spatially resolved observations of these 
disks, much of our knowledge of their structure is still de- 
rived from spatially Mnresolved spectroscopy and spectral 
energy distributions (SEDs). The interpretation of this in- 
formation (as well as spatially resolved data) requires the 
use of theoretical models, preferentially with as much re- 
alism and self-consistency as possible. Such disk models 
have been developed and improved over many years. When 
they are in reasonable agreement with observations they can 
also serve as a background onto which other processes are 
modeled, such as chemistry, grain growth, and ultimately 
the formation of planets. 

This chapter reviews the development of such self- 
consistent disk structure models, and discusses the cuiTent 
status of the field. We focus on the regions of the disk 
beyond the dust sublimation radius, since the very inner 
regions are discussed in the chapter by Najita et al.. To 
limit our scope further, we restrict our review to models 
primarily aimed at a comparison with observations. We 



will start with a concise resume of the formation and vis- 
cous evolution of disks (Section 2). This sets the radial disk 
structure as a function of time. We then turn our attention 
to the vertical structure, under the simplifying assumption 
that the gas temperature equals the dust temperature every- 
where (Section 3). While this assumption is valid in the 
main body of the disk, it breaks down in the disk surface 
layers. Section 4 treats the gas physics and chemistry of 
these surface layers, where much of the spectra originate. 
Photoevaporation flows also originate from the warm sur- 
face layers, and affect the disk evolution, which is the topic 
of Section 5. 

2. FORMATION AND VISCOUS EVOLUTION OF 
DISKS 

The formation of stars and planetary systems starts with 
the gravitational collapse of a dense molecular cloud core. 
Since such a core will always have some angular momen- 
tum at the onset of collapse, most of the infalling matter 
will not fall directly onto the protostar, but form a disk 
around it (e.g., Terebey et al, 1984; Yorke et al., 1993) 
or fragment into a multiple stellar system (e.g., Matsumoto 
and Hanawa, 2003). Because of the complexity of the lat- 
ter, we focus on the single star formation scenario here. 
While matter falls onto the disk, viscous stresses and grav- 
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Fig. 1. — Evolution of various disk and star quantities as a 
function of time after the onset of collapse of the cloud core 
(after Hueso and Guillot, 2005). SoUd lines: stellar mass 
(upper) and disk mass (lower). Dotted line: accretion rate 
in the disk. In this model the disk is formed aXt~ 0.03 
Myr, causing the jump in the dotted line at this point. The 
collapse phase is finished by 2 x 10'' years. 

itational torques within the disk wiU transport angular mo- 
mentum to its outer regions. As a consequence of this, 
most of the disk matter moves inward, adding matter to the 
protostar, while some disk matter moves outward, absorb- 
ing all the angular momentum (Lynden-Bell and Pringle, 
1974). During its formation and evolution a disk will spread 
out to several 100 AU or more (Nakamoto and Nakagawa, 
1994, henceforth NN94; Hueso and Guillot, 2005, hence- 
forth HG05). This spreading is only stopped when pro- 
cesses such as photoevaporation (this chapter), stellar en- 
counters (Scally and Clarke, 2001; Pfalzner et ai, 2005) or 
a binary companion (Artymowicz and Lubow, 1994) trun- 
cate the disk from the outside. During the collapse phase, 
which lasts a fewx 10"'' years, the accretion rate within the 
disk is very high {M ^ 10"'"^ - 10"'' Mq/yv), but quickly 
drops to M ^ 10"^ - 10~^ Mq/yt once the infall phase 
is over (NN94, HG05). The optical and ultraviolet excess 
observed from classical T Tauri stars (CTTSs) and Herbig 
Ae/Be stars (HAeBes) confirms that this on-going accre- 
tion indeed takes place {Calvet et al, 2000, and references 
therein). In Fig. 1 we show the evolution of various disk 
and star parameters. 

2.1. Anomalous viscosity 

An issue that is still a matter of debate is what consti- 
tutes the viscosity required for disk accretion, particularly 
after infall has ceased. Molecular viscosity is too small 
to account for the observed mass accretion rates. Turbu- 
lent and magnetic stresses, however, can constitute some 
kind of anomalous viscosity. The magnetorotational insta- 
bility (MRI), present in weakly magnetized disks, is the 
most accepted mechanism to drive turbulence in disks and 
transport angular momentum outwards (Balbus and Haw- 
ley, 1991; Stone and Pringle, 2001; Wardle, 2004, and ref- 



erences therein). 

There is a disk region (0.2 < r < 4 AU, for typi- 
cal CTTS disk parameters according to D'Alessio et al, 
1998) in which the ionization fraction is smaller than the 
minimum value required for the MRI. Neither thermal ion- 
ization (requiring a temperature higher than 1000 K), cos- 
mic ray ionization (requiring a mass surface density smaller 
than ~ 100 g/cm^) (Jin, 1996; Gammie, 1996), nor X-rays 
(Glassgold et al., 1997a, 1997b) are able to provide a suffi- 
cient number of free electrons to have MRI operating near 
the midplane. Gammie (1996) proposed a layered accretion 
disk model, in which a "dead zone" is encased between two 
actively accreting layers. The precise extent of this dead 
zone is difficult to assess, because the number density of 
free electrons depends on detailed chemistry as well as the 
dust grain size distribution, since dust grains tend to capture 
free electrons (Sano et al., 2000, and references therein). If 
the disk dust is like the interstellar dust, the MRI should 
be inhibited in large parts of the disk (Ilgner and Nelson, 
2006), though this is still under debate (e.g. Semenov et al, 
2004; see chapter by Bergin et al.). 

There are also other (non-magnetic) mechanisms for 
anomalous viscosity, like the baroclinic instability (Klahr 
and Bodenheimer, 2003) or the shear instability (Dubrulle 
et al, 2005), which are still subject to some controversy 
(see the recent review by Gammie and Johnson, 2005). An- 
gular momentum can also be transferred by global torques, 
such as through gravitational spiral waves (Tohline and 
Hachisu, 1990; Laughlin and Bodenheimer, 1994; Pickett 
et al, 2003 and references therein) or via global magnetic 
fields threading the disk (Stehle and Spruit, 2001), possibly 
with hydromagnetic winds launched along them (Blandford 
and Payne, 1982; Reyes-Ruiz and Stepinski, 1996). 

2.2. a-Dlsk models for protoplanetary disks 

To avoid having to solve the problem of viscosity in 
detail, but still be able to produce sensible disk mod- 
els, Shakura and Sunyaev (1973) introduced the "a- 
prescription", based on dimensional arguments. In this 
recipe the vertically averaged viscosity u at radius r is writ- 
ten as ly = aHpCs, where Hp is the pressure scale height 
of the disk and Cg is the isothermal sound speed, both eval- 
uated at the disk midplane where most of the mass is con- 
centrated. The parameter a summarizes the uncertainties 
related to the sources of anomalous viscosity, and is often 
taken to be of the order of a ~ 10~^ for sufficiently ionized 
disks. 

From conservation of angular momentum, the mass sur- 
face density S of a steady disk (i.e. with a constant mass 
accretion rate M), for radii much larger than the disk in- 
ner radius ri„, can be written as S(r) ~ M/Sttv. With 
Hp = Cs/flK, where CIk is the Keplerian angular velocity, 
we see that for r rin, 



where T'c(r) is the midplane temperature of the disk 
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at radius r and K is a constant with the value K = 
\/GM^fimp/3iTk. Here fi is the mean molecular weight in 
units of the proton mass rup, G is the gravitational constant, 
k is Boltzmann's constant and M^, is the stellar mass. As 
we will show in Section 3, most of the disk is 'irradiation- 
dominated', and consequently has temperature given ap- 
proximately by Tc ^ r~^/^. This results in the surface 
density going as E ~ r~^. This surface density distribu- 
tion is less steep than the so-called "minimum mass solar 
nebula" (MMSN), given by E ~ r~^/^ (Weidenschilling, 
1977; Hayashi, 1981). Strictly, the MMSN does not neces- 
sarily represent the mass distribution at any instant, but the 
minimum mass that has passed through the disk during its 
lifetime {Lissauer, 1993, and references therein). 

In reahty protoplanetary disks are not quite steady. After 
the main infall phase is over, the disk is not supplied any- 
more with new matter, and the continuing accretion onto the 
star will drain matter from the disk (see Fig. 1). In addition 
the disk viscously expands and is subject to photoevapora- 
tion (see Section 5). The timescale for "viscous evolution" 
depends on radius and is given by fyis — /v, which for 
typical CTTS parameters is 1 Myr at r ~ 100 AU. Since 
for irradiated disks tvis r, the outer regions evolve the 
slowest, yet they contain most of the mass. These regions 
(> 50 — 100 AU) therefore form a reservoir of mass con- 
stantly resupplying the inner regions. The latter can thus be 
approximately described by steady accretion disk models. 

There might also be dramatic variabiUty taking place 
on shorter timescales, as shown by FU Ori and EX Lupi 
type outbursts (Gammie and Johnson, 2005, and refer- 
ences therein). These outbursts can have various trigger- 
ing mechanisms, such as thermal instability (Kawazoe and 
Mineshige, 1993; Bell and Lin, 1994); close passage of 
a companion star {Bonnell and Bastien 1992; Clarke and 
Syer, 1996); mass accumulation in the dead zone followed 
by gravitational instability (Gammie, 1996; Armitage et al., 
2001). Disks are therefore quite time- varying, and constant 
a steady disk models should be taken as zeroth-order esti- 
mates of the disk structure. 

Given the challenges of understanding the disk viscos- 
ity from first principles, attempts have been made to find 
observational constraints on disk evolution {Ruden and Pol- 
lack, 1991; Cassen, 1996; Hartmann et al., 1998; Stepinski, 
1998). For example, Hartmann et al. (1998) study a large 
sample of CTTSs and find a dechne in mass accretion rate 
with time, roughly described as M ~ t~^'^, which they 
compare to the analytic similarity solutions of Lynden-Bell 
and Pringle (1974) for the expanding disk. A similar type 
of observational constraint is the recently found rough cor- 
relation M cx Ml (MuzeroUe et al, 2003a, 2005; Natta 
et al., 2004). High angular resolution mm- wave continuum 
imaging can also help to constrain the mass surface density 
distribution. With this technique Wilner et al. (2000) con- 
cluded that E oc for TW Hydra. Kitamura et al. (2002) 
find E r~P, with p = — 1 for a sample of T Tauri stars. 



3. VERTICAL STRUCTURE OF DUSTY DISKS 

With the radial structure following from accretion physics, 
as described above, the next issue is the vertical structure 
of these disks. Many authors have modeled this with full 
time-dependent 2D/3D (magneto/radiation-) hydrodynam- 
ics (e.g.. Boss, 1996, 1997; Yorke and Bodenheimer, 1999; 
Fromang et al., 2004). While this approach is obviously 
very powerful, it suffers from large computational costs, 
and often requires strong simplifying assumptions in the ra- 
diative transfer to keep the problem tractable. For compari- 
son to observed spectra and images these models are there- 
fore less practical. Most observation-oriented disk structure 
models split the disk into a series of (nearly independent) 
annuli, each constituting a 1-D or a two-layer local vertical 
structure problem. In this section we review this kind of 
'1+lD' models, and their 2-D/3-D generalizations. 

3.1. Basic principles 

The main objective of the models described in this section is 
the determination of the density and temperature structure 
of the disk. For a given surface density E(r), and a given 
gas temperature structure Tg{r^ z) (where z is the vertical 
coordinate measured upward from the midplane) the verti- 
cal density distribution p(r, z) can be readily obtained by 
integrating the vertical equation of hydrostatics: 

dP 

^ = -P^'k^ (2) 

where P = pc^ with = kTg/ jinip. The main com- 
plexity of a disk model Ues in the computation of the tem- 
perature structure. Since the main source of opacity is the 
dust, most models so far make the assumption that the gas 
temperature is equal to the dust temperature, so that the gas 
temperature determination reduces to solving a dust contin- 
uum radiative transfer problem. In Section 4 we will relax 
this assumption, but until then we will keep it. 

The temperature of the disk is set by a balance between 
heating and coohng. The disk cools by thermal emission 
from the dust grains at infrared wavelengths. This radia- 
tion is what is observed as infrared dust continuum radiation 
from such disks. Line cooling is only a minor coolant, and 
only plays a role for Tg when gas and dust are thermally 
decoupled. Dust grains can be heated in part by radiation 
from other grains in the disk. The iterative absorption and 
re-emission of infrared radiation by dust grains in the disk 
causes the radiation to propagate through the disk in a diffu- 
sive way. Net energy input comes from absorption of direct 
stellar Ught in the disk's surface layers, and from viscous 
dissipation of gravitational energy in the disk due to accre- 
tion. For most disks around CTTSs and Herbig Ae/Be stars 
the heating by stellar radiation is dominant over the viscous 
heating (except in the very inner regions). Only for strongly 
accreting disks does the latter dominate. 

Once the temperature structure is determined, the SED 
can be computed. The observable thermal emission of a 
dusty disk model consists of three wavelength regions (see 
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Fig. 2. — Build-up of the SED of a flaring protoplanetary 
disk and the origin of various components: the near-infrared 
bump comes from the inner rim, the infrared dust features 
from the warm surface layer, and the underlying continuum 
from the deeper (cooler) disk regions. Typically the near- 
and mid-infrared emission comes from small radii, while 
the far-infrared comes from the outer disk regions. The 
(sub-)millimeter emission mostly comes from the midplane 
of the outer disk. Scattering is not included here. 

Fig. 2). The main portion of the energy is emitted in a wave- 
length range depending on the minimum and maximum 
temperature of the dust in the disk. We call this the "ener- 
getic domain" of the SED, which typically ranges from 1 .5 
^m to about 100 /im. At shorter wavelength the SED turns 
over into the "Wien domain". At longer wavelengths the 
SED turns over into the "Rayleigh-Jeans domain", a steep, 
nearly powerlaw profile with a slope depending on grain 
properties and disk optical depth (see chapter by Natta et 
al.). Differences in disk geometry are mainly reflected in 
the energetic domain of the SED, while the submm and mm 
fluxes probe the disk mass. 

3.2. A first confrontation witli observations 

It is quite challenging to solve the entire disk structure ac- 
cording to the above principles. Early disk models were 
therefore often based on strong simplifications. An exam- 
ple of such a model is a perfectly flat disk being irradiated 
by the star due to the star's non-negligible size (Adams and 
Shu, 192,6; Fried] ung, 1985). The stellar radiation impinges 
onto the flat disk under an irradiation angle y ~ 0.4r*/r 
(with the stellar radius). Neglecting viscous dissipation, 
the effective temperature of the disk is set by a balance be- 
tween the irradiated flux (l/2)(^L*/47rr^ (withi* the stel- 
lar luminosity) and blackbody cooling crT^g, which yields 
Teff (X r~^/^. The energetic domain of its SED therefore 
has a slope of uFn oc with s = 4/3 = 1.33, which 



foUows from the fact that any disk with Tes oc has 
an SED slope of s = {4q — 2)/q. This steep slope arises 
because most of the steUar radiation is absorbed and re- 
emitted at small radii where the disk is hot. This produces 
strong emission at short wavelength. The long wavelength 
flux is weak because only little stellar radiation is absorbed 
at large radii. Observations of CTTSs, however, show SED 
slopes typically in the range s = 0.6 to 1 (Kenyan and Hart- 
mann, 1995), i.e. much less steep. The SEDs of Herbig 
Ae/Be stars show a similar picture, but with a somewhat 
larger spread in s, though it must be kept in mind that the 
determination of the slope of a bumpy SED hke in Fig. 2 is 
somewhat subjective. Meeus et al. (2001, henceforth MOl) 
divide the SEDs of Herbig Ae/Be stars into two groups: 
those with strong far-infrared flux (called 'group I', hav- 
ing slope s ~ —1 ... 0.2) and those with weak far-infrared 
flux (called 'group IF, having slope s ^ 0.2 ... 1). All 
but the most extreme group II sources have a slope that is 
clearly inconsistent with that of a flat disk. Note, at this 
point, that the Meeus 'group I' and 'group IF are unrelated 
to the Lada 'class F and 'class IF classification (both Meeus 
group I and II are members of Lada class II). 

A number of authors have employed another model to in- 
terpret their observations of protoplanetary disks: that of a 
steady accretion disk heated by viscous dissipation (Rucin- 
ski et al., 1985; Bertout et al., 1988; Hillenbrand et al., 
1992). These models are based on the model by Shakura 
and Sunyaev (1973). A detailed vertical structure model of 
such a disk was presented by Bell et al. (1997). The lumi- 
nosity of such disks, including the magnetospheric accre- 
tion column, is iaccr = GMt,M /vt:. For r » rin the 
effective temperature of such disks is given by crT^fg = 
3Mr2|./87r (with a the Stefan-Boltzmann constant), yield- 
ing an SED slope of s = 4/3, like for passive flat disks 
(Lynden-Bell, 1969; see solid lines of Fig. 6). Therefore 
these models are not very succesful either, except for mod- 
eling very active disks like FU Orionis (FUor) outbursts 
(see Bell and Lin, 1994). 

3.3. Flaring disk geometry 

It was recognized by Kenyan and Hartmann (1987) that a 
natural explanation for the strong far-infrared flux (i.e. shal- 
low SED slope) of most sources is a flaring ("bowl-shaped") 
geometry of the disk's surface, as depicted in Fig. 2. The 
flaring geometry allows the disk to capture a significant por- 
tion of the stellar radiation at large radii where the disk is 
cool, thereby boosting the mid- to far-infrared emission. 

The flaring geometry adds an extra term to the irradiation 
angle: ~ 0.4 Tt^jr^ rd{Hslr)/dr (Chiang and Goldre- 
ich, 1997, henceforth CG97), where Hs is the height above 
the midplane where the disk becomes optically thick to the 
impinging stellar radiation. In the same way as for the flat 
disks the thermal balance determines the Tea of the disk, but 
this now depends strongly on the shape of the disk: Hgir). 
The pressure scale height Hp, on the other hand, depends on 
the midplane temperature Tc by Hp = ^/ kTcV^ / jimpGM^ 
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(with the stellar mass). If we set Tc = T^s and if the 
ratio x = Hs/ Hp is known, then the system of equations is 
closed and can be solved (see appendix Chiang etai, 200 1 ) . 
For the special case that x is constant we obtain Hg oc r^/^, 
a-posteriori confirming that the disk indeed has a "bowl" 
shape. In general, though, x must be computed numeri- 
cally, and depends on the dust opacity of the disk upper 
layers. The resulting temperature profile is typically about 

The total luminosity of such a non-accreling flaring disk 
is Ldisk = C L^,, where C is the covering fraction of the 
disk. The covering fraction is the fraction of the starUght 
that is captured by the material in the disk. With a large 
enough disk inner radius (rjn 3> r*) one can write C ~ 
ma.x(Hs{r)/r). For an infinitely thin disk extending from 
r, to r — > 00 one has C = 0.25. The observed flux ratio 
Fdisk/F^ may deviate from C by about a factor of 2 due to 
the anisotropy of disk emission. 

In addition to irradiation by the star, the outer regions of 
a strongly accreting flaring disk (like an FUor object) can 
also be irradiated by the accretion luminosity from the inner 
disk {Kenyon and Hartmann, 1991; Bell, 1999; Lachaume, 
2004) and by the emission from the magnetospheric accre- 
tion column or boundary layer (Muzerolle et al, 2003b). 

3.4. Warm dust surface layer 

A closer look at the physics of an irradiation-dominated 
disk (be it flat or flared) reveals that its surface tempera- 
ture is generally higher than its interior temperature (Calvet 
etai, 1991; Malbet and Bertout, 1991; CG97). Dust grains 
in the surface layers are directly exposed to the stellar ra- 
diation, and are therefore hotter than dust grains residing 
deep in the disk which only 'see' the infrared emission by 
other dust grains. The temperature difference is typically 
a factor of 2 - 4 for non/weakly-accreting disks (see curve 
labeled "-9" in Fig. 4). For non-negligible accretion, on the 
other hand, the disk is heated from inside as well, producing 
a temperature minimum somewhere between the equatorial 
plane and the surface layer (see other curves in Fig. 4). Be- 
cause of the shallow incidence angle of the stellar radia- 
tion (f <^1, the vertical optical depth of this warm surface 
layer is very low. The layer produces optically thin emis- 
sion at a temperature higher than the effective temperature 
of the disk. The frequency-integrated flux of this emission 
is the same as that from the disk interior. As a consequence, 
the thermal radiation from these surface layers produces 
dust features in emission. This is exactly what is seen in 
nearly all non-edge-on T Tauri and Herbig Ae/Be star spec- 
tra (e.g. MOl; Kessler-Silacci et al, 2006), indicating that 
these disks are nearly always dominated by irradiation. 

3.5. Detailed models for flaring disks 

5.5.7. Disk structure. Armed with the concepts of disk 
flaring and hot surface layers, a number of authors pub- 
lished detailed 1+lD disk models and two-layer (sur- 
face+interior) models with direct appUcability to observa- 
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Fig. 3. — Characteristic radii of a disk around a 0.9 Mq star 
with T* = 4000 K and i?* = 1.9 at different accretion 
rates. The silicate sublimation radius is the dust 'irmer rim' . 
Figure based on models \yj D'Alessio et al. (1998). 

tions. The aforementioned CG97 model (with refinements 
described in Chiang et al., 2001) is a handy two-layer 
model for the interpretation of SEDs and dust emission 
features from non-accreting ('passive') disks. Lachaume et 
al. (2003) extended it to include viscous dissipation. 

The models by D'Alessio et al. (1998) solve the com- 
plete 1+lD disk structore problem with diffusive radiative 
transfer, including stellar irradiation and viscous dissipation 
(using the a prescription). The main input parameters are a 
global (constant) mass accretion rate M and a. The surface 
density profile E(r) is calculated self-consistently. This 
model shows that the disk can be divided into three zones: 
an outer zone in which the disk is dominated by irradiation, 
an inner zone where viscous dissipation dominates the en- 
ergy balance, and an intermediate zone where the midplane 
temperature is dominated by the viscous dissipation but the 
surface temperature by irradiation (See Fig. 3). In the inter- 
mediate zone the vertical disk thickness is set by M and a 
but the infrared spectrum is stiU powered by irradiation. In 
Fig. 4 the vertical structure of the disk is shown, for fixed S 
but varying M for r — \ AU. 

The models described by Dullemond, et al. (2002) apply 
exact 1-D wavelength-dependent radiative transfer for the 
vertical structure, but these models do not include viscous 
dissipation. 

3.5.2. Dust growth and sedimentation. Models of the kind 
discussed above describe the SEDs of CTTSs reasonably 
well. However, D'Alessio et al. (1999) argue that they tend 
to slightly overproduce far-infrared flux and have too thick 
dark lanes in images of edge-on disks. They also show 
that the percentage of expected edge-on disks appears to be 
overpredicted. They suggest that dust sedimentation could 
help to solve this problem. Chiang et al. (2001) find sim- 
ilar results for a subset of their Herbig Ae/Be star sample: 
the Meeus group 11 sources (see also CG97). They fit these 
sources by mimicking dust settling through a reduction of 
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Fig. 4. — Vertical temperature distribution of an irradiated 
a-disk at 1 AU, for a fixed S (chosen to be that of a disk 
model with M = 10~^ M© /yr for a = 0.01), but varying 
M, computed using the models of D'Alessio et al. (1998). 
The labels of the curves denote the 10-log of the accretion 
rate in Mq/yi. 

the disk surface height. Self-consistent computations of 
dust sedimentation produce similar SEDs and confirm the 
dust settling idea (Miyake and Nakagawa, 1995; Dullemond 
and Dominik, 2004b, henceforth DD04b; D'Alessio et al, 
2006). The disk thickness and far-infrared flux can also be 
reduced by grain growth {D 'Alessio et al. , 200 1 ; Dullemond 
and Dominik, 2004a). The chapter by Dominik et al. dis- 
cusses such models of grain growth and sedimentation in 
detail. 

From comparing infrared and (sub-)millimeter spectra of 
the same sources {Acke et al, 2004), it is clear that smaU 
and big grains co-exist in these disks. The (sub-)milUmeter 
spectral slopes usually require mm-sized grains near the 
midplane in the outer regions of the disk, while infrared 
dust emission features clearly prove that the disk surface 
layers are dominated by grains no larger than a few mi- 
crons (see chapter by Natta et al.). It appears that a bi- 
modal grain size distribution can fit the observed spectra: 
a portion of sub-micron grains in the surface layers respon- 
sible for the infrared dust enoission features and a portion 
of mm-sized grains in the disk interior accounting for the 
(sub-)miUimeter emission {Natta et al., 2001). 

3.6. The dust 'inner rim' 

The very inner part of the disk is dust-free due to dust 
sublimation (see chapter by Najita et al. for a discussion of 
this region). The dusty part of the disk can therefore be ex- 
pected to have a relatively abrupt inner edge at about 0.5 
AU for a 50 L© star (scaling roughly with \/TJl). If the gas 
inward of this dust inner rim is optically thin, which appears 
to be mostly the case (Muzerolle etal, 2004), then this dust 
iimer rim is illuminated by the star at a ~ 90 degree angle, 
and is hence expected to be much hotter than the rest of 
the disk behind it which is irradiated under a shaUow angle 



<C 1. {Natta et al, 2001). Consequently it must be hy- 
drostatically 'puffed-up', although this is still under debate. 
Natta et al. (2001) showed that the emission from such a 
hot inner rim can explain the near-infrared bump seen in 
almost all Herbig Ae/Be star SEDs (see e.g., MOl). This 
is a natural explanation, since dust sublimation occurs typ- 
ically around 1500 K, and a 1500 K blackbody bump fits 
reasonably well to the near-infrared bumps in those sources. 
Tuthill et al. (2001) independently discovered a bright half- 
moon ring around the Herbig Be star LkHa-101, which they 
attribute to a bright inner disk rim due to dust sublima- 
tion. Dullemond et al. (2001; henceforth DDNOl) extended 
the CG97 model to include such a puffed-up rim, and Do- 
minik et al. (2003) showed that the Meeus sample of Her- 
big Ae/Be stars can be reasonably well fitted by this model. 
However, for Meeus group II sources these fits required rel- 
atively small disks (see, however. Section 3.7). 

The initial rim models were rather simphfied, treating 
it as a vertical blackbody 'wall' (DDNOl). Isella and 
Natta (2005) improved this by noting that the weak de- 
pendence of the sublimation temperature on gas density is 
enough to strongly round off the rim. Rounded-off rims ap- 
pear to be more consistent with observations than the verti- 
cal ones: their flux is less inclination dependent, and their 
images on the sky are not so much one-sided. There is stiU 
a worry, though, whether the rims can be high enough to flt 
sources with a strong near-infrared bump. 

With near-infrared interferometry the rim can be spa- 
tially resolved, and thus the models can be tested. The 
measurements so far do not yet give images, but the mea- 
sured 'visibilities' can be compared to models. In this way 
one can measure the radius of the rim (e.g., Monnier et al., 
2005; Akeson et al, 2005) and its inclination (e.g., Eisner 
et al, 2003). Moreover it can test whether indeed the near- 
infrared emission comes from the inner rim of the dust disk 
in the first place (some doubts have been voiced by Vinkovic 
et al., 2003). We refer to the chapter by MiUan-Gabet et al. 
for a more in-depth discussion of interferometric measure- 
ments of disks. 

The inner rim model has so far been mainly appUed to 
Herbig Ae/Be stars because the rim appears so apparent 
in the near-infrared (NIR). But Muzerolle et al. (2003b) 
showed that it also applies to T Tauri stars. In that case, 
however, the luminosity from the magnetospheric accre- 
tion shock is required in addition to the stellar luminosity 
to power the inner rim emission. 

In addition to being a strong source of NIR flux, the 
'puffed-up' inner dust rim might also be responsible for 
the irregular few-day-long extinction events observed to- 
ward UX Orionis stars {Natta et al. 2001; Dullemond et 
al. 2003). The latter authors argued that this only works for 
self-shadowed (or only weakly flaring) disks (see Section 
3.7). 

In Fig. 5 we summarize in a qualitative way how the disk 
geometry (inner rim, flaring) affects the SED shape of an ir- 
radiated passive disk. In Fig. 6 the SEDs of actively accret- 
ing disks are shown, in which the irradiation by the central 
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Fig. 5. — Overall SED shape for non-accreting disks with 
stellar irradiation, computed using the 2-D radiative trans- 
fer tools from Dullemond & Dominik (2004a). The stellar 
spectrum is added in grey-scale. Scattered light is not in- 
cluded in these SEDs. SoUd line is normal flaring disk with 
inner dust rim; dashed line is when the rim is made higher; 
dot-dashed line is when the flaring is reduced (or when the 
disk becomes 'self-shadowed'); dotted line is when the in- 
ner rim is at 10 x larger radius. 



star is ignored. In reality, both the accretional heating and 
the irradiation by the central star must be included in the 
models simultaneously. 

3.7. 2-D radiative transfer in disk models 

The models described so far are all based on an approxi- 
mate 1h-1D (or two-layer) irradiation-angle description. In 
reality the structure of these disks is 2-D, if axisymmetry 
can be assumed, and 3-D if it cannot. Over the last 10 years 
many multi-dimensional dust continuum radiative transfer 
programs and algorithms were developed for this purpose 
(e.g., Whitney et al, 1992; Lucy et al, 1999; Wolf et al, 
1999; Bjorkman and Wood, 2001; Nicolinni et al, 2003; 
Steinacker et al, 2003). Most applications of these codes 
assume a given density distribution and compute spectra 
and images. There is a vast literature on such applications 
which we will not review here (see chapter by Watson et 
al.). But there is a trend to include the self-consistent ver- 
tical density structure into the models by iterating between 
radiative transfer and the vertical pressure balance equation 
(Nomura, 2002; Dullemond, 2002, henceforth D02; Dulle- 
mond and Dominik, 2004a, henceforth DD04a; Walker et 
al., 2004). The main improvements of 2-D/3-D models over 
1h-1D models is their ability to account for radial radiative 
energy diffusion in the disk, for cooling of the outer disk in 
radial direction, for the complex 3-D structure of the dust 
inner rim, and in general for more realistic model images. 

In addition to this, 2-D/3-D models allow for a 'new' 
class of disk geometries to be investigated. The 1-i-lD mod- 
els can, because of their reliance on an irradiation angle 
ip, only model disk geometries that are either flat or flared. 



Fig. 6. — Overall SED shape for accreting disks with- 
out stellar irradiation for two accretion rates. A sim- 
ple Shakura-Sunyaev model is used here with grey opac- 
ities. Solid line: pure Shakura-Sunyaev model (star not 
included); dotted hue: model with disk-self-irradiation in- 
cluded; dashed line: model with disk- self-irradiation and 
irradiation by the magnetospheric accretion column on the 
star included. 

In principle, however, there might be circumstances under 
which, roughly speaking, the surface of the outer disk re- 
gions Ues within the shadow of the inner disk regions (al- 
though the concept of 'shadow' must be used with care 
here). These shadowed regions are cooler than they would 
be if the disk was flaring, but the 2-D/3-D nature of radia- 
tive transfer prevents them from becoming entirely cold and 
flat. For Herbig Ae/Be stars the origin of the shadow might 
be the puffed-up inner rim (D02, DD04a), while for T Tauri 
stars it might be the entire inner flaring disk region out to 
some radius (DD04b). 

Although the concept of 'self- shadowing' is still un- 
der debate, it might be linked to various observable fea- 
tures of protoplanetary disks. For instance, DD04a showed 
that self-shadowed disks produce SEDs consistent with 
Meeus group II sources, while flaring disks generally pro- 
duce group I type SEDs, unless the disk outer radius is 
very small. It might also underly the observed correla- 
tion between SED shape and sub-millimeter slope (Acke et 
al. 2004). Moreover, self-shadowed disks, when spatially 
resolved in scattered light, would be much dimmer than flar- 
ing disks. 

4. GAS TEMPERATURE AND LINE SPECTRA 

Although the dust in disks is generally more easily ob- 
served, there is an obvious interest in direct observations 
of the gas. It dominates the mass, sets the structure and im- 
pacts dust dynamics and settling in these disks. Moreover, it 
is important to estimate how long disks remain gas-rich, and 
whether this is consistent with the formation time scale of 
gas giant planets (Hubickyj et al, 2004). Unfortunately, gas 
lines such as CO rotational, H2 rotational and atomic fine 
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structure lines often probe those surface regions of disks in 
which the gas temperature is difficuh to compute. The disk 
models we described above assume that the gas temperature 
in the disk is always equal to the local dust temperature. 
While this is presumably true for most of the matter deep 
within optically thick disks, in the tenuous surface layers of 
these disks (or throughout optically thin disks) the densities 
become so low that the gas will thermally decouple from 
the dust. The gas will acquire its own temperature, which 
is set by a balance between various heating- and cooUng 
processes. These processes depend strongly on the abun- 
dance of various atomic and molecular species, which, for 
their part, depend strongly on the temperature. The gas tem- 
perature, density, chemistry, radiative transfer and radiation 
environment are therefore intimately intertwined and have 
to be studied as a single entity. This greatly complicates 
the modeling effort, and the first models which study this in 
detail have only recently been published. 

This chapter focuses on stationary models, i.e. models 
that are in chemical, thermal and hydrostatic equilibrium. 
For the tenuous regions of disks the chemical time scales 
are short enough that this is valid, in contrast to the longer 
chemical time scales deeper in the disk (e.g., Aikawa and 
Herbst, 1999; Willacy etal, 2000). The models constructed 
so far either solve the gas temperature/chemistry for a 
fixed gas density structure (Jonkheid et ai, 2004; Kamp 
and Dullemond, 2004), or include the gas density in the 
computation to obtain a self-consistent thermo-chemical- 
hydrostatic structure (Gorti and Hollenbach, 2004; Nomura 
and Millar, 2005). 

4.1. Basic gas physics 

The physics and chemistry of the surface layers of pro- 
toplanetary disks strongly resembles that of photon dom- 
inated regions (PDRs, Tielens and Hollenbach, 1985; Ya- 
mashita et al. 1993). In those surface layers the gas temper- 
ature generally greatly exceeds the dust temperature. But 
the dust-gas coupling gradually takes over the gas tempera- 
ture balance as one gets deeper into the disk, typically be- 
yond a vertical column depth of Ay — 1, and forces the gas 
temperature to the dust temperature. 

The uppermost surface layer contains mostly atomic and 
ionized species, since the high UV irradiation effectively 
dissociates all molecules (Aikawa et al., 2002). The photo- 
chemistry is driven by the stellar UV irradiation and/or in 
case of nearby O/B stars, by external illumination. In flar- 
ing disk models, the stellar UV radiation penetrates the disk 
under an irradiation angle ^p hke the one described in the 
previous section. This radiation gets diluted with increas- 
ing distance from the central star and attenuated by dust 
and gas along an inclined path into the disk. The stellar UV 
radiation therefore penetrates less deep into the disk than 
external UV radiation. As one goes deeper into the surface 
layer, the gas becomes molecular (see chapter by Bergin et 
al.). 

The thermal balance of the gas in disks is solved by 



equating all relevant heating and cooling processes. For this 
gas thermal balance equation, a limited set of key atomic 
and molecular species is sufficient: e.g., H2, CO, OH, H2O, 
C+, O, Si+ and various other heavy elements. For most 
atoms and molecules, the statistical equilibrium equation 
has to include the pumping of the fine structure and ro- 
tational levels by the cosmic background radiation, which 
become important deep in the disk, where stellar radiation 
cannot penetrate. The full radiative transfer in chemical 
models is very challenging, and therefore generally approx- 
imated by a simple escape probability approach, where the 
pumping and escape probability are derived from the opti- 
cal depth of the line (similar to the approach of Tielens and 
Hollenbach, 1985 for PDRs). Even though the emitted pho- 
tons travel in all directions, the optical depth used for this 
escape probability is the line optical depth in the vertical 
direction where the photons most readily escape. 

One of the most critical ingredients of these models is the 
UV and X-ray radiation field (stellar and external), which 
can be split into the far-ultraviolet (FUV, 6-13.6 eV), the 
extreme-ultraviolet (EUV, 13.6-100 eV) and X-ray (> 100 
eV) regime. In the literature the far ultraviolet radiation 
field (FUV) is often represented by a single parameter Go 
describing the integrated intensity between 912 and 2000 A 
normalized to that of the typical interstellar radiation field 
(Habing, 1968). However, several papers have shown the 
importance of a more detailed description of the radiation 
field for calculations of the chemistry and the gas heat- 
ing/cooling balance (Spaans et al., 1994; Kamp & Bertoldi, 
2000; Bergin et al., 2003; Kamp et al, 2006; Nomura and 
Millar, 2005). For instance, in T Tauri stars the radiation 
field is dominated by strong Ly a emission, which has con- 
sequences for the photodissociation rate of molecules that 
can be dissociated by Ly a photons. The photoelectric heat- 
ing process, on the other hand, depends strongly on the 
overall shape of the radiation field, which is much steeper 
in the case of cool stars. A similar problem appears in the 
X-ray spectra of cool M stars, which are dominated by Une 
emission. 

FUV induced grain photoelectric heating is often a dom- 
inant heating process for the gas in the irradiated surface 
layers. The FUV photon is absorbed by a dust grain or a 
polycyclic aromatic hydrocarbon (PAH) molecule, which 
ejects an energetic electron to the gas, and heats the gas via 
the thermalization of the energetic electron. Its efficiency 
and thus the final gas temperature depends strongly on the 
grain charge, dust grain size and composition (PAHs, sil- 
icates, graphites, ices, etc.). X-rays from the central star 
also heat only the uppermost surface layers, as their heating 
drops off monotonically with column, and gets quite small 
by columns of order 10^^ cm^^. 

4.2. Surfaces of optically thick disks 

This subsection focuses on the warm surface layers of 
the optically thick disk at Ay < 1, measured verticaUy 
downwards, where gas and dust temperatures decouple. 
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Fig. 7. — Upper panel: gas (solid Une) and dust (dashed 
line) temperatures in a vertical cut through the T Tauri 
star model at 100 AU. Overplotted are the most important 
chemical transitions from atomic to molecular species. The 
black filled circle shows the point where a photoevapora- 
tive flow can be initiated (see Section 5). Lower panel: 
same vertical cut for the hydrogen number density (fixed 
input density distribution). Overplotted are the ranges over 
which the respective heating/cooling processes contribute 
more than 10% to the total heating/cooUng rate at that depth 
('PE' means photoelectric heating; 'H2,torm' means heating 
through molecular hydrogen formation; see Section 5 for 
the definition of Tevap.crit and TviriaO- 



Modeling of these surface layers is not affected by the op- 
tically thick disk interior and depends mainly on the local 
UV/X-ray flux and the gas density. The typical hydrogen 
gas number density in these layers is roughly nniAy = 
1) ~ 10^ (100 AU/r) S-^ cm-3. The location of the 
Ay = 1 surface depends on the ratio 5 of the dust surface 
area per hydrogen nucleus to the interstellar value, which 
is roughly 10~^^cm^/H. Assuming a surface density that 
drops linearly with radius and E(l AU) 1000 g cm^^, 
the fractional column density Sgurf/S contained in the sur- 
face layer (Tgas ^ Tdust) is usually smaU, Ssurf/S ~ 
1.5 X 10~^ (r/AU), but increases linearly with radius. 

4.2.1. Gas temperatures. The detailed temperature struc- 
ture of the surface layers of optically thick young disks was 
studied for the first time by Jonkheid et al. (2004), Kamp 
and Dullemond (2004), and Nomura and Millar (2005). 
These models neglect EUV irradiation and start with neu- 
tral hydrogen in the top layers. Fig. 7 shows the vertical 



structure in a disk model with 0.01 Mq at 100 AU around 
a 0.5 Mq T Tauri star (from models of Kamp and Dulle- 
mond, 2004; note that the density structure in these models 
is not iterated with the gas temperature). Very high in the 
atmosphere at particle densities as low as n < 10'' cm~^ 
(Ay < 10~^), the gas temperature is set by a balance be- 
tween photoelectric heating and fine structure line cooUng 
of neutral oxygen (Kamp and Dullemond, 2004; Jonkheid 
et al, 2004). This leads to gas temperatures of several hun- 
dred K. Deeper in the disk, for Ay > 0.01, molecules can 
shield themselves from the dissociating FUV radiation. As 
soon as the fraction of molecular hydrogen becomes larger 
than 1%, H2 Une cooUng becomes important. Molecular 
line emission - mainly CO and H2 - cools the gas down 
to below hundred K before the densities are high enough 
for gas and dust to thermally couple. As gas temperatures 
drop below ^ 100 K, H2 no longer contributes to the cool- 
ing. Instead CO, which has a rich rotational spectrum at 
low temperatures, becomes an important coolant. At larger 
radii the FUV flux from the central star drops as well as the 
density of the surface layer, leading to lower gas tempera- 
tures. At distances r > 100 AU the gas temperature is too 
low for the endothermic destruction of H2 by O atoms and 
hence the surface layer at those distances contains substan- 
tial fractions of molecular hydrogen. 

4.2.2. Implications for the disk structure. Detailed mod- 
els of the gas temperature have shown that gas and dust are 
collisionally coupled at optical depth Ay > 1. Thus the 
basic assumption Tgas = Tkust of the disk structure mod- 
els presented in the previous section is justified for the disk 
interior The main effect of the higher gas temperatures in 
the warm surface layer is an enhanced flaring of the disk 
{Nomura and Millar, 2005). 

4.2.3. Observations and comparison with models. The 
pure rotational lines of H2 such as J = 2 - S(0) [28 /im], J = 
3-1 S(l)[17/xm],J = 4-2S(2)[12^m]andJ = 6-4S(4) 
[8 /im] trace the warm gas (100-200 K) in the disks. Even 
though there is some controversy about detection of those 
lines with different instruments (Thi et al., 2001; Richter 
et al, 2002), there is a tentative detection of H2 in AB Au- 
rigae using the Texas Echelon Cross Echelle Spectrograph 
(TEXES) at the Infrared Telescope Facility {Richter et al., 
2002). Bary et al. (2003) report v = 1 - S(l) |2.12 /im] 
emission in high resolution spectra {R ^ 60 000) of the T 
Tauri stars GG Tau A, LkCa 15, TW Hya and DoAr21 . This 
emission most Ukely arises in the low density, high temper- 
ature upper surfaces beyond 10 AU. According to the disk 
models, warm H2 exists indeed in the opticaUy thin surface 
layers, where Tgas ^ ^dust and the observed fluxes can 
be reproduced {Nomura and Millar, 2005). Fig. 8 reveals 
the effect of UV fluorescence on the line strength. The UV 
fluorescent lines, which are an excellent probe of the inner 
disk (r < few AU), are discussed in detail in the chapter 
by Najita et al.. The detection of the mid-lR H2 lines at 
low spectral resolution (e.g., with Spitzer) is hindered by 
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Fig. 9. — Impact of detailed gas temperature modeling 
and dust settling on the emission lines of oxygen and CO 
{Jonkheid et at. 2004): Tgas = Tdust (dotted line), detailed 
gas energy balance (solid line), dust settling (dashed line). 



the low Une-to-continuum ratio. 

The impact of detailed gas modeling differs for the var- 
ious emission lines. CO, which forms deeper in the disk is 
generally less affected than fine structure lines such as [O i] 
and [C II] that form in the uppermost surface layers, where 
Tgas 3> Tdust (Fig. 9). Since the gas temperature in those 
layers is set by photoelectric heating, dust settling leads to 
lower temperatures and thus to weaker line emission. 

The [O l] 6300 A line is another tracer of the physics in 
the tenuous surface layers (see also the chapter by Bergin 
et al.). It has been detected in a number of externally il- 
luminated proplyds in the Orion nebula (Johnstone et al, 
1998) as well as in T Tauri and Herbig Ae/Be stars (Acke 
et al., 2005; Acke and van den Ancker, 2006). Storzer and 
Hollenbach (1999) explain the emission in the Orion pro- 
plyds by the photodissociation of the OH molecule, which 
leaves about 50% of the atomic oxygen formed in the up- 



per level of the 6300 A line. Acke et al. (2005) find 
indication of Keplerian rotation from the [O l] line profiles. 
However, they need OH abundances higher than those pre- 
dicted from disk models to fit the emission from the disks 
around Herbig Ae/Be stars. Gas models of those disks re- 
veal the presence of a high temperature reservoir (few 1000 
K); hence the [O i] line might arise partly from thermal ex- 
citation at radii smaller than 100 AU (Kamp et al., 2006). 
Resolved [Oi] 6300 A line emission from the disk around 
the Herbig Ae star HD 100546 (Acke and van den Ancker, 
2006) shows that the emission is spread between ^ 1 and 
100 AU and supports the presence of a gap at ~ 10 AU as 
reported initially by Bouwman et al. (2003). 

4.3. Optically thin disks 

As protoplanetary disks evolve, the dust grains grow to at 
least centimeter sizes and the disks become optically thin. 
In addition, as we shall discuss Section 5, the gas in the disk 
ultimately disappears, turning the disk into a debris disk. It 
is therefore theoretically conceivable that there exists a tran- 
sition period in which the disk has become optically thin 
in dust continuum, but still contains a detectable amount of 
gas. The source HD141569 (5 Myr) might be an example of 
this, as Brittain et al. (2003) observed UV excited warm CO 
gas from the inner rim at ~ 17 AU, and Dent et al. (2005) 
cold gas further out (J = 3 - 2). Measuring the gas mass 
in such transition disks sets a timescale for the planet for- 
mation process. The Spitzer Legacy Science Program 'For- 
mation and Evolution of Planetary Systems' (FEPS) has set 
upper Umits on gas masses of ~ 0.1 Mj around solar-type 
stars with ages greater than 10 Myr (Meyer et al., 2004; 
Hollenbach et al., 2005; Pascucci et al., 2005). 

4.3.1. Disk models. Several groups have so far studied 
these transition phases of protoplanetary disks: Gorti and 
Hollenbach (2004) modeled the disk structure and gas/dust 
emission from intermediate aged disks around low-mass 
stars, Kamp and Bertoldi (2000), Kamp and van Zadel- 
hoff (2001), and Kamp et al. (2003) modeled the gas chem- 
istry and line emission from A-type stars such as /3 Pictoris 
and Vega. Jonkheid et al. (2006) studied the gas chemi- 
cal structure and molecular emission in the transition phase 
disk around HD 141569 A. These models are all based on 
the same physics as outlined above for the optically thick 
protoplanetary disks. The disks are still in hydrostatic equi- 
librium, so that the disk structure in these low mass disks is 
similar to that in the more massive disks with the midplane 
simply removed. However, some fundamental differences 
remain: the minimum grain size in these disks is typically 
a few microns, much larger than in the young protoplane- 
tary disks; in addition, the dust may have settled towards 
the midplane, and much of the solid mass may reside in 
larger particles (a. > 1 cm) than can be currently observed. 
This reduces the grain opacity and the dust-to-gas mass ra- 
tio compared to the younger optically thick disks. At ra- 
dial midplane gas column densities smaller than lO^'' cm^^, 
these disks are optically thin to stellar UV and ~ 1 keV X- 
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Fig. 10. — Mid-infrared spectrum in the 24 - 40 /xm wave- 
length region showing the dust continuum and dominant gas 
emission lines for a disk model with Mgas = 10~^ Afj and 
-^dust = 10~^ Mj (with dust defined as particles smaller 
than 1 nrni). A distance to the disk of 30 pc and a spectral 
resolving power R = 600 is assumed {Gorti and Hollen- 
bach, 2004) 

ray photons and the gas is mostly atomic. At radial columns 
greater than that, the gas opacity becomes large enough 
to shield H2 and CO, allowing significant molecular abun- 
dances. For disks extended to 100 AU, very Uttle mass (very 
roughly > 10~^ Mj) is needed to provide this shielding. 

4.3.2. Comparison with observations. Since the contin- 
uum remains optically thin, the mid-lR spectrum is dom- 
inated by fine-structure emission lines from ions such as 
[Fe II] and [Si 11]; large columns of neutral sulphur are com- 
mon, leading to strong [S l] emission (Fig. 10). However, 
the strength and thus detectability of these lines depends on 
the abundances of heavy metals in late phases of disk evolu- 
tion, which is uncertain, especially the more refractory Fe. 
In the opaque molecular regions somewhat closer to the star, 
the gas temperature exceeds 100 K and molecular hydrogen 
emission is produced. The S(0) and S(l) H2 lines stay opti- 
cally thin over a large range of disk masses and if detected 
are more diagnostic of disk mass than other fine structure 
lines {Gorti and Hollenhach (2004)). While the strongest 
molecular bands of H2O (important coolant in disk mid- 
plane) and OH are similar in strength to the fine structure 
lines, the H2 lines are weak and detection can be signifi- 
cantly hampered by the low line-to-continuum ratio (weak 
narrow line against the bright dust thermal background). 
These mid-IR lines generally originate from 1-10 AU. 

Kamp et al. (2003) have shown that beyond 40 AU the 
dominant coolant for the latest tenuous stages of disk evo- 



lution is the [Cii] 158 iim Une. The fine structure lines 
of C, O and C+ trace only the surface of these tenuous 
disks: [Oi] becomes rapidly optically thick and C"*" and 
C turn into CO as soon as UV CO and H2 bands become 
optically thick, and stellar UV cannot penetrate any further. 
Since typical gas temperatures are higher than in molecular 
clouds, CO lines from the upper rotational levels (J = 4 - 
3) are predicted to be stronger than the lower J lines. Dent 
et al. (2005) have recently detected the CO J = 3 - 2 line 
in HD141569 and disk modeling by Jonkheid et al. (2006) 
shows that the profile excludes a significant contribution 
from gas inwards of ~ 80 AU and estimate the total gas 
mass to be 80 Mg. 

5. PHOTOEVAPORATION OF A DISK BY ITS 
CENTRAL STAR 

5.1. Introduction 

The above section has shown that in the surface layers of 

the disk the gas temperature can become very high, greatly 
exceeding the dust temperature. The warm surface gas can 
flow off the disk and escape the gravity of the star. Since 
the heating process responsible for these high temperatures 
is the radiation from the central star or a nearby O-star, this 
process is caUed "photoevaporation". The viscous evolu- 
tion (i.e. accretion and spreading) of the disk, discussed in 
Section 2, can be strongly affected by this photoevaporation 
process. TypicaUy, it significantly shortens the 'Ufetime' of 
a disk compared to pure viscous evolution. Photoevapo- 
ration can also create inner holes or truncate the outer disk. 
This has relevance to observations of such disks, such as the 
percentage of young stars with infrared excess versus their 
age (Haisch et al, 2001; Carpenter et al., 2005), or the in- 
ferred 'large inner holes' of some disks (e.g., Calvet et al., 
2002; Bouwman et al., 2003; Forrest et al., 2004; D'Alessio 
et al., 2005). It has also far-reaching consequences for the 
formation of planets, as we wiU discuss below. 

Photoevaporation has already been discussed in earlier 
reviews (Hollenbach et al., 2000; Hollenbach and Adams, 
2004; Richling etal., 2006). However, these reviews mainly 
focused on the heating by a nearby massive star (such as the 
famous case of the proplyds in Orion). In contrast, in this 
section we will exclusively review recent results on pho- 
toevaporation by the central star, consistent with the pre- 
vious sections which focus on heating and photodissocia- 
tion by the central star Progress in this field since PPIV 
has been mostly theoretical, since observations of diagnos- 
tic gas spectral lines for the case of photoevaporation by the 
central, low mass star requires greater sensitivity, spectral 
resolution, and spatial resolution than currently available. 
We will, however, discuss the implications for the observed 
'irmer holes' and disk lifetimes. 

5.2. The Physics of Photoevaporation 

5.2.1. Basic Concepts. Photoevaporation results when 

stellar radiation heats the disk surface and resulting thermal 
pressure gradients drive an expanding hydrodynamical flow 
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to space. As shown in Section 4 the main heating photons 
he in the FUV, EUV and X-ray energy regimes. X-rays, 
however, were shown to be of lesser importance for pho- 
toevaporation {Alexander et al, 2004b), and we will not 
consider them further. 

There are two main sources of the strong EUV and FUV 
excesses observed in young low mass stars: accretion lumi- 
nosity and prolonged excess chromospheric activity. Recent 
work {Alexander et al. , 2004a) has shown that EUV photons 
do not penetrate accretion columns, so that accretion cannot 
provide escaping EUV photons to power photoevaporation. 
Alexander et al. (2005) present indirect observational ev- 
idence that an active chromosphere may persist in T Tauri 
stars even without strong accretion, and that EUV luminosi- 
ties of $EUV > 10"^^ photons/s may persist in low mass 
stars for extended (> 10^'^ — lO*^ yrs) periods to illuminate 
their outer disks. FUV photons may penetrate accretion 
columns and also are produced in active chromospheres. 
They are measured in nearby, young, solar mass stars with 
little accretion and typically (with great scatter) have lumi- 
nosity ratios I/FuvZ-^'boi ~ 10~^ or $fuv ~ 10**^ pho- 
tons/s. 

EUV photons ionize the hydrogen in the very upper lay- 
ers of the disk and heat it to a temperature of ~ 10** K, inde- 
pendent of radius. FUV photons penetrate deeper into the 
disk and heat the gas to T ^ 100 — 5000 K, depending on 
the intensity of the FUV flux, the gas density and the chem- 
istry (as was discussed in Section 4). Whether the EUV 
or FUV heating is enough to drive an evaporative flow de- 
pends on how the resulting thermal speed (or sound speed) 
compares to the local escape speed from the gravitationally 
bound system. A characteristic radius for thermal evapora- 
tion is the "gravitational radius" Vg, where the sound speed 
equals the escape speed: 

^« = T^~'°°^^ ilOOOKj iMoj-^'^ 

Early analytic models made the simple assumption that 
photoevaporation occurred for r > rg, and that the warm 
surface was gravitationally bound for r < rg. However, 
a closer look at the gas dynamics shows that this division 
happens not at rg but at about 0.1 - 0.2 {Lijfman, 2003; 
Adams et al., 2004; Font et al, 2004), and that this division 
is not entirely sharp. In other words, photoevaporation hap- 
pens mostly outside of the "critical radius" r„ ~ O.lSr^, 
though a weak evaporation occurs inside of rcr- Since these 
are important new insights since PPIV, we devote a subsub- 
section on them below. 

Witii T ~ 10^ K tiie critical radius for EUV-induced 
photoevaporation is (EUV) 1 - 2{M^/Mq ) AU. 
However, there is no fixed rcr (FUV) because the FUV- 
heated gas has temperatures that depend on FUV flux and 
gas density, i.e., on r and z. Therefore, rcr(FUV) depends 
on r and z, and may range from 3-150 AU for solar mass 
stars. 

The evaporative mass flux S depends not only on the 
temperature of the photon-heated gas, but also on the ver- 



tical penetration depth of the FUV/EUV photons. For 
EUV photons this is roughly set for r < rcr ^ 1 AU 
by the Stromgren condition that recombinations in the 
ionized layer equal the incident ionizing flux. Neglect- 
ing dust attenuation, this penetration column can be ex- 
pressed: Ay (EUV) - 0.05 5$4(^(r/AU)-i/^ where 
$41 = ^Euv/W^ photons/s and S is the ratio of the dust 
surface area per hydrogen to the interstellar dust value (see 
Section 4 and note that S can be much smaller than unity if 
dust has settled or coagulated). Outside of 1 AU, the pene- 
tration depth falls even faster with r, roughly as r^"^/^ (see 
Hollenbach et al. 1994). On the other hand, the FUV pene- 
tration depth is set by dust attenuation, or Ay (FUV) ^ (p, 
where we recall that (p is the irradiation angle and depends 
on disk flaring. In general Ay (EUV) < Ay (FUV), so flie 
EUV-ionized skin of the disk lies on top of the FUV-heated 
gas surface layer. 

The penetration depth is an important quantity because 
it sets the density at the base of the photoevaporative flow: 
the deeper the penetration depth, the higher the density. The 
flux of outflowing matter is proportional to the product of 
local density and sound speed within this heated layer. This 
is why the complex surface structure models of Section 
4 are so important for FUV-driven photoevaporation. For 
EUV-driven photoevaporation, on the other hand, the situa- 
tion is less complicated, since the temperature in the ionized 
skin of the disk is independent of r and z, as long as 2; > Zb, 
where Zb is the bottom of the ionized layer, i.e. the base of 
the flow. For this simple case, the evaporative mass flux 
originates at Zb, which is where the highest density gas at 
temperature Teuv — 10* K resides. 

Although FUV-heated layers have lower temperatures 
than the EUV-heated skin they are at higher densities and 
may equally well initiate the flow and determine the mass 
flux as EUV photons (see Johnstone et al, 1998 for a sim- 
ilar situation for externally illuminated disks). Gorti and 
Hollenbach (in preparation, henceforth GH06) find that the 
FUV-photoevaporative flow typically originates at vertical 
heights where T ~ 100 - 200 K, yielding rcr ~ 50 - 100 
AU. For r > 50 AU, the FUV photoevaporation dominates. 
On the other hand, EUV photons (with rcr ^ 1 AU) affect 
the planet forming regions at r <C 50 AU more than the 
FUV photons. 

5.2.2. Photoevaporation as a Bernoulli flow. One way 
to understand why the disk can evaporate at radii as small 
as 0.2rg is to consider the evaporative flow as a Bernoulli 
flow along streamlines {Lijfman, 2003; Adams et al., 2004). 
These streamlines initially rise nearly vertically out of the 
disk and then bend over to become asymptotically radially 
outward streamhnes. If a streamline starts at r > rg, then 
the flow rapidly goes through a sonic point and achieves the 
sound speed Cg near the base of the flow. The mass flux rate 
in the flow is then S ~ pbCs, where pb is the mass density 
of the gas at the base. 

On the other hand, if a streamline starts at r <C rg, 
the gas at its base lies deep in the gravitational potential. 
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As a simplification let us now treat these streamlines as 
if they are entirely radial streamlines (ignoring their ver- 
tical rise out of a disk). Then the standard atmospheric 
solution has a density that falls off from r to roughly 
as exp(— rg/2r). The gas flows subsonically and acceler- 
ates, as it slowly expands outward, until it passes through 
a sonic point at rg < 0.5rg (O.brg is the classic Parker 
wind solution for zero rotation). For r <^ Vg, the mass 
flux is reduced considerably by the rapid faU-off of the den- 
sity from r to r^. For r < rg, the mass flux is roughly 
given by the density at times the sound speed times the 
dilution factor (rg/r)^ that accounts for mass conservation 
between r and Tg: S ~ pbe~^^^'^^Cs{rs/r)'^. Assuming the 
same pb and Cg at all r, we see that E(0.2rg) ^ 0.5S(rg) 
and that S(O.lrg) ~ 0.17E(rg). This demonstrates that 
Tcr ~ 0.1 5rg for this simpUfied case, and that even for 
r < evaporation is weak, but not zero. In Fig. 7 the base 
of the flow is marked with the large dot (though that figure 
shows a static, non-evaporating model with only FUV heat- 
ing). In that figure, Tviriai is the temperature such that the 
sound speed equals the escape speed; Tevap.crit = 0.2Tviiiai 
is roughly where the photoevaporation flow originates (i.e., 
where r = rcr). 

5.2.3. Mass loss rates for EUV-induced flows. Although 
central star FUV models are not yet published, several cen- 
tral star EUV models have appeared in the hterature. Hol- 
lenbach et al. (1994) first outhned the essential physics of 
EUV-induced flows by the central star and presented an ap- 
proximate analytic solution to the mass loss rate for a disk 
larger than rg. The basic physics is the Stromgren relation, 
$EUV — ar'nf.r'^, where is the hydrogen recombination 
coefficient and rig is the electron density in the ionized sur- 
face gas. This sets the hydrogen nucleus (proton) number 

1/2 

density at the base of the flow rib 'I'euV' ^'^'^ therefore 
an identical proportionality for the mass loss rate: 

*-~4xl0-."(^)°'(^)"^ (4, 

in units of Mq /yr. Radiation hydrodynamical simulations 
(Yorke and Welz, 1996; Richling and Yorke, 1997) find a 
similar power-law index for the dependence of the mass- 
loss rate on the EUV photon rate of the central star. This 
result applies for both high and low mass central stars, and 
is valid for a weak stellar wind. The effect of a strong stellar 
wind is such that the ram pressure reduces the scale height 
of the atmosphere above the disk and the EUV photons are 
allowed to penetrate more easily to larger radii. This in- 
creases the mass-loss rate from the outer parts of the disk. 
It is noteworthy that the diffuse EUV field, caused by re- 
combining electrons and protons in the disk's ionized at- 
mosphere inside rcr, controls the EUV-induced mass-loss 
rates (Hollenbach et al., 1994) for disks with no or small 
iimer holes (< rcr). This effect negates any potential for 
self-shadowing of the EUV by the disk. 



5.3. Evolution of photoevaporating disks 

5.3.1. Case without viscous evolution. Let us first assume 
a disk that does not viscously evolve: it just passively un- 
dergoes photoevaporation. For disks with size ra < rcr, the 
photoevaporation proceeds from outside in. The mass flux 
rate at r^ is much higher than inside of r^, because the gas 
at rd is least bound. In addition, disk surface densities gen- 
erally fall with r (see Section 2). Therefore, the disk shrinks 
as it photoevaporates, and most of the mass flux comes from 
the outer disk radius. However, for disks with r^ > rcr, two 
types of disk evolution may occur. For EUV photons, Hol- 
lenbach et al. (1994) showed that the mass flux S beyond 
rcr goes roughly as r~'^^^ if there is no inner hole extending 
to rcr- The timescale for complete evaporation at r goes as 
S(r)/E(r). As long as S does not drop faster than r^^"'', 
the disk will evaporate first at r ~ rcr, and, once a gap 
forms there, will then steadily erode the disk from this gap 
outwards. 

If, on the other hand, I](r)/i](r) decreases with r, then 
the disk shrinks from outside in as in the rj, < r^ case. 
The photoevaporation by the FUV from the central star has 
not yet been fully explored, but preliminary work by GH06 
suggests that the mass flux S in the outer disks around solar 
mass stars increases with r. In this case, the disk evaporates 
from outside in for most generally assumed surface density 
laws, which decrease with r. The combined effect of EUV 
and FUV photoevaporation then is likely to erode the disk 
outwards from rcr (EUV) ~ 1 AU by the EUV flow and 
inwards from the outer disk radius by the FUV flow, sand- 
wiching the intermediate radii. 

5.3.2. Case with viscous evolution. Now let us consider a 
disk that is actively accreting onto the star (see Section 2). 
In general, if the photoevaporation drills a hole somewhere 
in the disk or 'eats' its way from outside in, the forces of 
viscous spreading tend to move matter toward these photo- 
evaporation regions, which can accelerate the dissipation 
of the disk. If the disk has a steady accretion rate M, 
then a gap forms once Mcvap oc r^S exceeds M. Since 
r^E oc r""-^ for EUV photoevaporation beyond rcr, the 
gap first forms at the minimum radius ( ) and then 
works its way outward. Clarke et al. (2001) presented time- 
dependent computations of the evolution of disks around 
low mass stars with $euv ~ 10^^^^'^ photons s^^. Their 
model combines EUV photoevaporation with a viscous evo- 
lution code. However, they used the current hypothesis at 
that time that evaporation only occurs outside of rg. After 
^ 10® to 10^ years of viscous evolution relatively unper- 
turbed by photoevaporation, the viscous accretion inflow 
rates fall below the photoevaporation rates at rg. At this 
point, a gap opens up at rg and the inner disk rapidly (on 
an inner disk viscous timescale of ^ 10'"' yr) drains onto 
the central star or spreads to rg where it evaporates. In this 
fashion, an inner hole is rapidly produced extending to rg . 

Alexander et al. (2006a, 2006b) extended the work of 
Clarke et al. to include the effect of rcr < Vg, and to treat 



13 



10^ 
10* 

102 

10° 

^ 10-2 
10-* 
10-6 



10-8 



: ^;5Myr 


Spreod'ng ond photoevoporotion ■ 
of T Touri disk- 










a 






6.06Myr | \ 6.12Myr " 


. . " 





0.1 



1.0 



10.0 
r [AU] 



100.0 



1000.0 



Fig. 11. — Evolution of the surface density of a EUV- 
photoevaporating disk (Figure adapted from Alexander et 
al., 2006b). This simulation starts from a given disk struc- 
ture of about 0.05 Mq (marked with 'Start' in the fig- 
ure). Initially the disk accretes and viscously spreads (soUd 
Unes). At t = 6 X 10^ yr the photoevaporation starts affect- 
ing the disk. Once the EUV-photoevaporation has drilled a 
gap in the disk at ~ 1 AU, the destruction of the disk goes 
very rapidly (dashed lines) .The inner disk quickly accretes 
onto the star, followed by a rapid erosion of the outer disk 
from inside out. In this model the disk viscosity spreads 
to > 1000 AU; however, FUV-photoevaporation (not in- 
cluded) will likely truncate the outer disk. 



the outward EUV evaporation of the disk beyond 1 
AU. They show that once the inner hole is produced, the dif- 
fuse flux from the atmosphere of the inner disk is removed 
and the attenuation of the direct flux by this same atmo- 
sphere is also removed. This enhances the EUV photoevap- 
oration rate by the direct EUV flux from the star, and the 

1 /2 

effect magnifies as the inner hole grows as AfEXJV « r-J^^^^, 
again derivable from a simple Stromgren criterion. The 
conclusion is that the outer disk is very rapidly cleared once 
the inner hole forms (see Fig. 1 1). 

The rapid formation of a cleared out inner hole almost 
instantly changes the nature and appearance of the disk. 
The above authors compare their model favorably with a 
number of observations: (i) the rapid transition from clas- 
sical T Tauri stars to weak line T Tauri stars, (ii) the al- 
most simultaneous loss of the outer disk (as detected by 
submillimeter measurements of the dust continuum) with 
the inner disk (as detected by near IR observations of very 
hot dust near the star), and (iii) the SED observations of 
large (3-10 AU) inner holes in those sources (see dotted 
line of Fig. 5) with evidence for low accretion rates and in- 
termediate mass outer disks such as the source CoKu Tau/4. 
Fig. 12 shows the evolutionary tracks of their models with 
^Euv = 10''^ photons/s compared to the observations of 
weak-line T Tauri stars (WTTSs) and CTTSs. 

In a similar vein Armitage et al. (2003) used the com- 
bination of EUV photoevaporation and viscous dispersal. 
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Fig. 12. — Near- to mid-infrared color (in magnitudes) 
versus 850 /im flux for photoevaporation/viscous evolution 
models. The data are taken from Hartmann et al. (2005) 
wA Andrews and Williams (2005): 850 /um detections (cir- 
cles) and upper limits (triangles) are plotted for both CTTSs 
(filled symbols) and WTTSs (open symbols). Evolution- 
ary tracks are shown for models with stellar masses 0.5 
(dashed), 1.0 (soHd), and 2.0 M0(dotted), at a disk incH- 
nation of i = 60° to the line of sight. The thick tracks 
to the right and left show the 1 Mq model at z = and 
i = 80°, respectively. Crosses are added every 1 Myr to 
show the temporal evolution. Initially the (optically thin) 
850 /Ltm flux declines slowly at constant (optically thick) 
infrared color However, once the viscous accretion rate 
falls below the photoevaporation rate, the disk is rapidly 
cleared from the inside-out. (Figure adapted from Alexan- 
der et al. 2006b.) 



together with an assumed dispersion of a factor of 3 in the 
initial disk mass, to explain the dispersion in the lifetime 
and accretion rates of T Tauri disks. They found that the 
models predict a low fraction of binaries that pair a classi- 
cal T Tauri star with a weak-lined T Tauri star. Their models 
are in better agreement with observations of disk lifetimes 
in binaries than models without photoevaporation. 

Going one step further, Takeuchi et al. (2005) con- 
structed models combining viscous evolution, EUV photo- 
evaporation, and the differential radial motion of grains and 
gas. Their models predicted that for low-mass stars with a 
low photoevaporation rate, dust-poor gas disks with an in- 
ner hole would form (WTTs), whereas for high mass stars 
(evolved Herbig Ae/Be) with a high photoevaporation rate, 
gas-poor dust rings would form. 

Matsuyama et al. (2003) pointed out that if the EUV lu- 
minosity is created by accretion onto the star, then, as the 
accretion rate diminishes, the EUV luminosity drops and 
the timescale to create a gap greatly increases. Even worse, 
as discussed above, the EUV photons are unlikely to es- 
cape the accretion column. Only if the EUV luminosity re- 
mains high due to chromospheric activity does EUV photo- 
evaporation play an important role in the evolution of disks 
around isolated low mass stars. Alexander et al. (2005) 



14 



argue this is the case. Ruden (2004) provides a detailed 
analytic analysis which describes the evolution of disks in 
the presence of viscous accretion and photoevaporation and 
compares his results favorably with these two groups. 

5.4. Effect on planet formation 

The processes which disperse the gas influence the for- 
mation of planets. Small dust particles follow the gas flow. 
If the gas is dispersed before the dust can grow, all the dust 
will be lost in the gas dispersal and planetesimals and plan- 
ets will not form. Even if there is time for particles to co- 
agulate and build sufficiently large rocky cores that can ac- 
crete gas {Pollack et al, 1996; Hubickyj et al, 2004), the 
formation of gas giant planets like Jupiter and Saturn will 
be suppressed if the gas is dispersed before the accretion 
can occur. Furthermore, gas dispersal helps trigger gravita- 
tional instabilities that may lead to planetesimal formation 
{Goldreich and Ward, 1973; Youdin and Shu, 2002; Throop 
etal, 2005), affects planet migration (e.g.. Ward, 1997) and 
influences the orbital parameters of planetesimals and plan- 
ets {Kominami and Ida, 2002). 

5.4.7. Gas Rich Versus Gas Poor Giant Planets in the So- 
lar System. Shu et al. (1993) showed that with $euv ~ 
10^^ photons s~^, the early Sun could have photoevapo- 
rated the gas beyond Saturn before the cores of Neptune 
and Uranus formed, leaving them gas poor. However, this 
model ignored photoevaporation inside of r^. The current 
work by Adams et al. (2004) would suggest rather rapid 
photoevaporation inside of 10 AU, and make the timing of 
this scenario less plausible. FUV photoevaporation (either 
from external sources or from the central star) may provide 
a better explanation. Prehminary results from GH06 sug- 
gest that the early Sun did not produce enough FUV gen- 
erally to rapidly remove the gas in the outer giant planet 
regions. Adams et al. and Hollenbach and Adams, (2005) 
discuss the external illumination case, which looks more 
plausible. 

5.4.2. Truncation of the Kuiper Belt. A number of ob- 
servations point to the truncation of Kuiper Belt Ob- 
jects (KBOs) beyond about 50 AU (e.g., Allen, Bernstein, 
and Malhotra, 2002; Trujillo and Brown, 2001). Adams 
et al. (2004) and Hollenbach and Adams (2004, 2005) show 
that photoevaporation by a nearby massive star could cause 
truncation of KBOs at about 100 AU, but probably not 50 
AU. The truncation is caused by the gas dispersal before 
the dust can coagulate to sizes which survive the gas dis- 
persal, and which can then later form KBOs. Models of 
FUV photoevaporation by the early Sun are needed. 

5.4.3. Formation of Planetesimals. In young disks, dust 
settles toward the midplane under the influence of the stellar 
gravity and coagulates. Once coagulated dust has concen- 
trated in the midplane, the roughly centimeter-sized parti- 
cles can grow further by collisions or by local gravitational 
instabiUty {Goldreich and Ward, 1973; Youdin and Shu, 



2002). A numerical model by Throop and Bally (2005) fol- 
lows the evolution of gas and dust independently and con- 
siders the effects of vertical sedimentation and external pho- 
toevaporation. The surface layer of the disk becomes dust- 
depleted which leads to dust-depleted evaporating flows. 
Because of the combined effects of the dust setthng and the 
gas evaporating, the dust-to-gas ratio in the disk midplane 
is so high that it meets the gravitational instability criteria 
of Youdin and Shu (2002), indicating that kilometer- sized 
planetesimals could spontaneously form. These results im- 
ply that photoevaporation may even trigger the formation of 
planetesimals. Presumably, photoevaporation by the central 
star may also produce this effect. 

6. SUMMARY AND OUTLOOK 

In this chapter we have given a brief outhne of how disks 
form and viscously evolve, what their structure is, what 
their spectra look like in dust continuum and in gas lines, 
and how they might come to their end by photoevaporation 
and viscous dispersion. The disk structure in dust and gas 
is summarized in Fig. 13. Evidently, due to the broadness 
of the topic we had to omit many important issues. For 
instance the formation of disks is presumably much more 
chaotic than the simple picture we have discussed. In re- 
cent years there is a trend to outfit even the observation- 
oriented workhorse models with ever more detailed physics. 
This is not a luxury, since the volume of observational 
data (both spectral and spatial) is increasing dramatically, 
as shown by various other chapters in this book. For in- 
stance, with observational information about dust growth 
and sedimentation in disks, it will be necessary to include 
realistic dust evolution models into the disk models. Ad- 
ditionally, with clear evidence for non-axial symmetry in 
many disks (e.g., Fukagawa et al, 2004) modelers may be 
forced to abandon the assumption of axial symmetry. The 
thermal-chemical study of the gas in the disk surface lay- 
ers is a rather new field, and more major developements are 
expected in the next few years, both in theory and in the 
comparison to observations. These new insights will also 
influence the models of FUV-photoevaporation, and thereby 
the expected disk lifetime. 

A crucial step to aim for in the future is the unification of 
the various aspects discussed here. They are all intimitely 
connected together and mutually influence each other. Such 
a unification opens up the perspective of connecting seem- 
ingly unrelated observations and thereby improving our un- 
derstanding of the bigger picture. 
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